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On the spectrum of a class of quantum models 
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PACS 03.65.Ge - Solutions of wave equations: bound states 
PACS 02 . 30 . Ik - Integrable systems 

PACS 42 . 50 . Pq - Cavity quantum electrodynamics; micromasers 

Abstract - The spectrum of any quantum model which eigenvalue equation reduces to a three- 
term recurrence, such as a displaced harmonic oscillator, the Jaynes-Cummings model, the Rabi 
model, and a generalized Rabi model, can be determined as zeros of a corresponding transcendental 
function F(x). The latter can be analytically determined as an infinite series defined solely in 
terms of the recurrence coefficients. 



Introduction. — The present work deals with a class 
1Z of quantum models described by a Hamiltonian H char- 
acterized in that the eigenvalue equation 



H4> = Ecj) 



(1) 



in the Bargmann Hilbcrt space of analytical functions B 
[UEj reduces to a three-term difference equation 



~t~ a n c n -\- b n c Ti 







(n>0). (2) 



Here {c n }^L are the sought expansion coefficients of a 
physical state described by an entire function 



n=0 



(3) 



The recurrence coefficients a n and b n are functions of 
model parameters and we require that b n =/= 0. In what 
follows, models of the class 1Z are assumed to satisfy 
lim„_i. 0O a n 7^ and limn-,.^ b n /a n = 0. The recurrence 
coefficients may have either finite limits 



a, 



(n 



(4) 



(a Poincare difference equation [3]), or a power-like depen- 
dence a n ~ an 4 , b n ~ bri" . 

The conventional boson annihilation and creation oper- 
ators d and a) satisfying commutation relation [a, af] = 1 
are represented in the Bargmann Hilbcrt space B as [HE] 

d 



a — > 



d- 



(5) 



Therefore, the class 1Z comprises among others single- 
mode boson models represented by 



H = Aitfa + A 2 a t + A 3 a + A 4 



(6) 



where Aj are in general matrix coefficients. Indeed, on 
comparing the same powers of z on both sides of ([T]), the 
consecutive by a, a', and d terms in lead to nc n , c n _i, 
and (n + l)c„+i terms in the recurrence (0, respectively. 
The eigenvalue equation ([T]) then reduces to a coupled sys- 
tem of ordinary differential equations [THUE] • The model 
Hamiltonian (|6]) encompasses many prominent examples, 
such as a displaced harmonic oscillator [T| , the Rabi model 
[T1l6ll7] . the Jaynes and Cummings (JC) model proposed as 
an approximation to the Rabi model [S] , and a generalized 
Rabi model [9] . The Rabi model describes the simplest in- 
teraction between a cavity mode with a bare frequency oj 
and a two-level system with a bare resonance frequency lo . 
The model is characterized by the Hamiltonian [UEIITHH] > 



hjja) a + \o~\(a) + a) + 



(7) 



where A is a coupling constant, and fi = Hljq/2. In what 
follows we assume the standard representation of the Pauli 
matrices <7j and set the Planck constant K = 1. 

Our main results is that the spectrum of any quantum 
model from 1Z can be obtained as zeros of a transcendental 
function defined by infinite series 



F(x) = a + 2J P1P2 ■■■Pk 



(8) 



fe=i 



solely in terms of the coefficients of the three-term recur- 
rence relation in eq. (El [10] : 



Pi 

u, 



pi = ui - 1, m = 1, 



I > 2. 



1 - ui-ibi/{aiai-{) 



(9) 
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Fig. 1: Fdho as a function of x = E/uj for « = \/uj = 0.7, 
A = jti/w = 0, in the case when the spectrum of the Rabi 
model reduces to that of a displaced harmonic oscillator. In 
agreement with the exact analytic formula (|15[) . the zeros of 
Fdho(x) satisfy xi = I — 0.49. 



The function F(x) is different from the G-functions of 
Braak [9]. The latter have only been obtained in the case 
of the Rabi model by making explicit use of the parity 
symmetry. In contrast to Braak's result [9], our F(x) can 
be straightforwardly defined for any model from 1Z. The 
definition of F(x) does not require either a discrete sym- 
metry or to solve the three-term difference equation @ 
explicitly. All what is needed to determine F(x) is an ex- 
plicit knowledge of the recurrence coefficients a n and b n . 
This brings about a straightforward numerical implemen- 
tation and results in a great simplification in determin- 
ing the spectrum. The proof of principle is demonstrated 
in fig. [1] in the case of exactly solvable displaced har- 
monic oscillator [T]. The ease in obtaining the spectrum 
is of importance regarding recent experimental advances 
in preparing (ultra) strongly interacting quantum systems 
[T^HTfj] , which can no longer be reliably described by the 
exactly solvable JC model [SJ, and wherein only the full 
quantum Rabi model can describe the observed physics. 

Proof of the main result. — In earlier studies fTJ 
[H[SJ|!], largely motivated by the Frobenius method [T7] 
of solving differential equations, the n = part of the 
recurrence ([2]) was taken as an initial condition. Given 
the requirement of analyticity [c-k = for k > 1; cf. cq. 
{3}], the three-term recurrence © degenerates for n = 
into an equation involving mere two terms and imposes 
that 

c 1 /c = -a . (10) 

By solving the recurrence upwardly, one as a rule arrives 
at the expansion coefficients of singular functions outside 
the physical Bargmann Hilbert space B [rHlfTT] . Only at a 
particular discrete set of parameter values that correspond 
to the spectrum of a model the functions would belong to 



B [TJ[S] . Our approach differs from earlier studies [T, 5 , 9 
in that our focus is initially on the truly three-term part 
of the recurrence ([2]) for n > 1. This would enable us to 
remain in the physical Bargmann Hilbert space B for any 
value of physical parameters. The remaining n — part 
of the recurrence © , which reduces to a relation between 
Ci and Co, then becomes the boundary condition for <fi G B 
defined by the n > 1 solutions of the recurrence ([2]). 

In virtue of the Perron and Kreuser generalizations 
(Theorems 2.2 and 2.3(a) in ref. [3]) of the Poincare the- 
orem (Theorem 2.1 in ref. [3J), the difference equation ^ 
(considered for n > 1) possesses in our case two linearly 
independent solutions: 

• (i) a dominant solution {dj}°^ and 

• (ii) a minimal solution {mj}j2. . 

The respective solutions differ in their limiting behavior 
c„+i/c„ in the limit n — > oo. The the so-called minimal 
solution guaranteed by Perron's and Kreuser's theorems 
satisfies [31IT5] 



lim m n+ i/m n = 0. 



(11) 



On substituting the minimal solution of @ for the c n 's in 
eq. (J3J), <j)(z) automatically becomes an entire function. In 
what follow, the entire function generated by the minimal 
solution of the n > 1 part of © will be denoted as (j) m (z). 
As demonstrated in the case of the Rabi model below, 
4> m (z) belongs to B irrespective of the value of the model 
parameters. 

Now it is important to realize that only the ratios of 
subsequent terms m n+ i/m n of the minimal solution (and 
hence not of the dominant solution) are related to infinite 
continued fractions (cf. Theorem 1.1 due to Pincherle in 
ref. [3J) 



m n+ i _ — b n+ i b n +2 b n+ 3 



))),, 



dn+l~ O-n+2— &n+3 — 



(12) 



Any minimal solution is, up to a multiplication by a con- 
stant, unique [3J. This has two immediate consequences. 
First, 4> m (z) is, up to a multiplication constant, unique. 
Second, the ratio ro = mi /too of the first two terms of a 
given minimal solution is unambiguously fixed. Note that 
the ratio ro = toi/too involves too, although it takes into 
account the recurrence ^ only for n > 1 [3J. The re- 
maining n = part of the recurrence ([2]) imposes another 
condition [cf. cq. (|T0|) ] on the ratio ro = mi/mo- The con- 
dition (|10p can be translated into the boundary condition 
on the logarithmic derivative of </> m (z), 



Mo) 



-a . 



(13) 



Obviously, one has ro =/= —ao in general. As a rule it is 
impossible to find a minimal solution satisfying an arbi- 
trarily prescribed initial condition on the ratio TOi/toq of 
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its first two terms. (This point was not made clear by ei- 
ther Schweber pQ or later on by Braak j5J.) If one defines 
the function F — oq + ro, with ro given by the continued 
fraction in ([T2"]) , the zeros of F would correspond to those 
points in a parameter space where the condition ([TU]) is 
satisfied. Then the coefficients of 4> m (z) £ B satisfy the 
recurrence @ including n = 0, and (f> m (z) belongs to the 
spectrum. Analogous to the Schrbdinger equation, the 
boundary condition (|13j) enforces quantization of energy 
levels. 

The rest of the proof follows upon an application of the 
Euler theorem, which transforms the continued fraction 
defining r in eq. (fl~2j) into an infinite series (cf. eqs. 
(4.4)- (4.5) on p. 43 of [3j forming the basis of the Hhird" 
method of Gautschi [3] at the special case of his N = 
that is also used in our numerical implementation [11)). 
The latter leads directly to cqs. © and © for F [TU] . 

Illustration of some basic properties. In the 

case of a displaced harmonic oscillator, which is the special 
case of Hr in ([7]) for \i = 0, the recurrence ([2]) becomes 
[cf. eq. (A.17) of ref. jT]] 



1 



Cn+l 



(n + 1)k 



1 



Cn-l = 0, 



(14) 



where dimcnsionless parameter k = A/a; reflects the cou- 
pling strength and x = e = E/lu is a dimensionless energy 
parameter. Energy levels satisfy [TJ 



ei = Ei/w = I - k 2 



(15) 



where I £ N is a nonnegative integer (including zero). As 
demonstrated in fig. [U Fdho{x) defined by eqs. ([8]) and 
([9]) displays for w = 1 and k = 0.7 a series of discontinuous 
branches extending monotonically between — oo and +oo, 

which intersect the x-axis at xi = I — 0.49, I = 0, 1, 

The zeros correspond exactly to the position of the energy 
levels (|T5|) . 

In the present example one can also explicitly illustrate 
that the difference equation with the condition (fTU|) 
imposed as an initial condition of the recurrence, and 
then solved upwardly, always possesses a solution. Obvi- 
ously, such a solution would be generally a dominant solu- 
tion. Consequently, the function <f>(z) defined by power se- 
ries expansion with the expansion coefficients {dj}°°^ 
would exhibit (typically branch-cut |T| ) singularities in the 
^-complex plane and would not belong to B. Indeed, the 
recurrence IT41) is solved for n > with 



i 4 a -*V), 



(16) 



where a = x + k 2 and L@ are associated Laguerre polyno- 
mials [TO] (note different sign of n compared to eq. (2.16) 
of Schweber [1]). The substitution (fT6)) transforms (jT4j) 
into a 3-point rule that is identically satisfied by the asso- 
ciated Laguerre polynomials |19j . The Rodrigues formula 

ESI 



i\ dz r ' 



(17) 



implies Lq'(z) = 1, L\ a ' (z) = —z + 1 + a, and in virtue 

of ma 



c = n a L[ a \n 2 ) = n c \ 

Cl = K a - 1 L[ a - 1 \ K 2 ) = K a ~ 1 x. 



(18) 



Thus ci/co = x/k, which is exactly the n = part of (I14|) . 
In view of the asymptotic [19] 



(n — > oo) 



where the (generalized) binomial coefficients 

a(a- l)(a-2)---(a-/c + l) 



one finds from (H 



Cn+l 



fc! 



(19) 



(20) 



la-n 



« n + 1 



(n -> oo) (21) 



whenever a ^ N. Therefore, unless a is a nonnegative 
integer, the solution of eq. (|16[) is the dominant solution, 
the power series solution ([3]) has only a finite radius of 
convergence and thus does not belong to B. Only if the 
condition f|15[) is satisfied, a = x + k 2 is an integer, the 
dominant solution changes smoothly into the minimal so- 
lution of the recurrence (|16[) . and the corresponding 
vanishes. 

In the case of the Rabi model, the recurrence © be- 
comes 

fn{x) 1 



Cn+l 



(n+1) 
where f n (x) is given by 



1 



0, 



/„(a:) = 2ft + 



1 

2ft 



(22) 



(23) 



k = A/a;, as in eq. (|T4^). and dimensionless A = fi/cj 
(cf. eq. (A8) of Schweber [T], which has mistyped sign 
in front of his 6„_i, and eqs. (4) and (5) of [9]). At a 
first glance the recurrence (|22j) does not reduce to (fT4]) for 
A = as one would expect from 0. The equivalence of 
eqs. (fT4]) and ([22]) for A = is disguised by the fact that 
the recurrence (j22[) has been obtained after the unitary 
transformation of Hr in ([7]) induced by D — exp[ft(a^ — 
S)<ti] JT] . Thereby the dimensionless energy parameter x 
in cqs. jUJ, ([23]) is a; = (E'/w) + ft 2 . 

One can verify that the entire function cj) m (z) generated 
by the minimal solution {m 3 -}^ belongs to the Bargmann 
Hilbcrt space of analytical functions B for all model pa- 
rameters, i.e. even if the quantization criterion F = is 
not satisfied. Linearly independent solutions of recurrence 
([221) can be written in the limit n > 1 as 



l/(n/c»), 
n x K n /Y(n 



dominant 
minimal. 



(24) 
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Let us consider the z — > oo asymptotic of 

oo 

X (z) = Y, v n(x)z n (25) 



n=0 



in the case of the minimal solution. Obviously, <j) m (z) ~ 
X(z) for z — > oo. The asymptotic of x( z ) can be deter- 
mined by the saddle point of the Euler-Maclaurin integral 
representation of x( z ) 



X {z)= u(s,x)e slnz ds + 0(\z\ a ) (z->oo), 



a being some constant, — 1 < a < and 



(26) 



y ( S)X ) M (— )V = e s [ ln(Ke )- lns ] +:rlns = e h{s) . (27) 

For z — > oo and | argz |< it/ 2 one finds a unique saddle 
point in the cut complex plane that completely governs 
the asymptotic behaviour of x( z ) 



X(z) 



-2tt 



Here 



h"(t s 



h"(t s ) 



1/2 



, K z+*[i+W"*)] + 0{\z\ a ). (28) 



1 



■ 7 ^2 ( 29 ) 

KZ + X (KZ + X) Z 

is the second derivative of h(s) defined by eq. ([27]) at the 
saddle point t s « kz + x. Hence x{ z ) behaves essentially 
as e KZ for z —¥ oo. Such a behaviour of x( z ) was largely to 
be expected in view of that the asymptotic of the minimal 
solution ([24)) is only a slightly perturbed version of 1 /T(n+ 

I)- 

Discrete symmetries. — The Rabi Hamiltonian Hr 
[eq. (|7J|] is known to possess a discrete Z2-symmctry cor- 
responding to the constant of motion, or parity, P = 
cxp(inJ) [HE], where 



a) a 



2 (! + 



(30) 



is the familiar operator known to generate a continuous 
U(l) symmetry of the JC model [HE]. In contrast to ref. 
[9], our approach does not necessitate any active use of 
any underlying discrete symmetry. This has been demon- 
strated in a recent comment |10j . where the regular spec- 
trum of the Rabi model has been reproduced as zeros of a 
corresponding Fnd(x) based on the recurrence ([2"2"]> . Ob- 
viously it was not possible then to determine what is the 
parity of a state corresponding to a given zero of Fr^x). 

Nevertheless, not only any discrete symmetry can be 
easily incorporated in our approach, this can be accom- 
plished more straightforwardly and more easy than in the 
approach of Braak (9j . Indeed, the Rabi Hamiltonian Hr 
is an example of a general Fulton and Gouterman Hamil- 
tonian of a two-level system [2"0] 



H FG = Al + Ba 1 +Ca 3 , 



(31) 



with 



A = luo) a, 



B = A(a t +a), C = fi. 



(32) 



The Fulton and Gouterman symmetry operation g [20) is 
realized by reflections 



a -+ —a, 



(33) 



which leave the boson number operator a) a invariant. Be- 
cause [g, A] = [g, C] = {g, B} = 0, go~ 3 is the symmetry of 
Hr [20) . To any cyclic Z2 operator, such as go's, one can 
associate a pair of projection operators 



P ± = \{i±go 3 ) 



(p ± ) 2 = p ± . 



(34) 



However, because of g<r 3 , the projectors P do not mix the 
upper and lower components of a wave function <fi (conven- 
tional Pauli representation of <j\ and o~ 3 is assumed). In 
order to employ the Fulton and Gouterman reduction in 
the positive and negative parity spaces, wherein one com- 
ponent of 4> is generated from the other by means of the 
operator g, one is forced to work in a unitary equivalent 
single-mode spin-boson picture 



oja'a + \io\ + A<73(a' + a). 



(35) 



The transformation is accomplished by means of the uni- 
tary operator U — (<7i + a 3 )/^/2 = U^ 1 . The transforma- 
tion interchanges the expressions for B and C in (|32[) . re- 
sulting in gu\ becoming the symmetry of Hr. Eqs. (4.12) 
and (4.13) of [20] then yield 



+ A(a f + a) ± ng]^ 



E 



(36) 



where the superscripts ± denote the positive and negative 
parity eigenstates of P^ [with a 3 being replaced by o\ in 
(|34p]. Working in the Bargmann space, 



[ujzd z + X(z + d z ) ± /i<?]0 d 



E 



(37) 



Eqs. ([37]) are equivalent to a coupled system of first-order 
eqs. (10) of the supplement to [9]. In contrast to [9], one 
does not need any ill motivated substitutions and vague 
arguments to arrive at the symmetry resolved spectrum 
of the Rabi model. Assuming power series expansions 
</> ± (z) = X^^Lo " 2 ™ f° r t ne positive and negative par- 
ity states, one arrives directly at the following three-term 
recurrence relations 



1 



-71+ 1 



n(n + 1) 



-[n-x±(-l) n A]<* 



1 



= 0, (38) 



where x = E^/ui and k = \/oj, as in eq. (j!4l) . and A = 
n/u>, as in eq. ([2"2"]l . In arriving at eqs. ([55]) we have 
merely used that 



--(z) = <p ± (-z) = j2(-iy 



c±z r - 



(39) 



71=0 
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Fig. 2: Fr(x) for k = 0.7, A = 0.4, and u = 1, i.e. the same 
parameters as in fig. 1 of ref. [9], shows zero at w —0.707805 
for a negative parity state and at w —0.4270437 for a positive 
parity state. After addition of k 2 = 0.49 they correspond to 
the zeros at —0.217805 and 0.0629563, respectively, in fig. 1 of 
refs. M\M- 



Upon defining corresponding F^(x) by means of eqs. ([5]) 
and ([9]) , one can not only recover the regular spectrum of 
the Rabi, but also distinguish different parity cigcnstatcs. 
This is demonstrated in fig. [2 

Fig. |3] demonstrates that on setting A = the parity 
eigenvalues become doubly-degenerate eigenvalues corre- 
sponding to those of a displaced harmonic oscillator. That 
was to be expected, because the three-term recurrence re- 
lations ([55)1 reduce to the three-term recurrence ([hi]) for 
A = 0. 

Discussion. — In his recent letter [9], Braak claimed 
to solve the Rabi model analytically (see also Viewpoint 
by Solano [24] ). He suggested that a regular spectrum of 
the Rabi model was given by the zeros of transcendental 
functions G±(x) in the variable x = (E/tu) + k 2 [9], 



G±(x) 



71 = 



K n (x) 



it 



A 



(40) 



Here the coefficients K n (x) were obtained recursively by 
solving the Poincare difference equation (|22|) upwardly 
starting from the initial condition (|10[) . Braak needed one 
page of arguments (between eq. (10) on p. 2 and the end 
of the first paragraph on p. 3 of the on-line supplement to 
[5]) to arrive at his G±{x). The arguments involved an ill 
motivated substitution and that a sufficient condition for 
the vanishing of an analytic function G+(x; z) (defined by 
cq. (16) of his supplement) for all z G C is if it vanishes 
at a single point z = 0. The argument is essential to ar- 
rive at (|40p . However, as an example of any homogeneous 
polynomial shows, the argument is obviously invalid. Ad- 
ditionally, the coefficients K n 's in (|4"0"|) at zeros of G±(x) 



Fig. 3: Energy levels for k — 0.7 and ui = 1 as a function of 
A evolve from the doubly degenerate level corresponding to 
those of the displaced harmonic oscillator at ti = / — 0.49 to 
the parity eigenstates of the Rabi model. 



have to be the minimal solutions. The difficulty of calcu- 
lating minimal solutions can be illustrated by an attempt 
to calculate the Bessel functions of the first kind J n (x) for 
fixed x = 1 by an upward three-term recurrence from the 
initial values Jo(x) and J\{x). It turns out that all digits 
calculated in single precision came out illusory already for 
n>7 (cf. Table 1 of [3]). 

We have advanced an alternative approach which does 
not require any explicit solution of any difference equation 
and which applies to an entire class of models. A symme- 
try is not required to characterize a quantum model in 
terms of a suitable F(x). Contrary to ref. [5], F(x) can 
be obtained directly from the recurrence coefficients char- 
acterizing an eigenvalue equation of a given model. In the 
case of the Rabi model, nothing but the elementary rela- 
tion (j39|) has been employed in arriving from (|37]l to ([38)1 , 
and the corresponding F^ (x) were unambiguously deter- 
mined from (|38[) through eqs. (J5J and ([9]). The resulting 
transcendental functions F^ (x) are different from G± (x) 
considered by Braak [9]. With the exception of the dis- 
placed harmonic oscillator and the JC model, the zeros of 
F(x) have to be still determined numerically. Therefore 
we disagree with Solano's view [24] that a mere characteri- 
zation of a model by a function such as G± (x) of ref. [9] , or 
our F(x), is tantamount to solving the model analytically 
in a closed-form. 

Our results rejuvenate the Schwcbcr quantization cri- 
terion ro + ao = known in the case when ro remained 
to be expressed in terms of continued fractions [e.g. eq. 
(fT2T) ] pQ. The Schweber criterion was deemed impractical 
and has not been employed to solve for the quantized en- 
ergy levels of any quantum model [UIHGOU] ■ The criterion 
was either regarded to require a numerical diagonalization 
in a truncated Hilbert space (pp. 4-5 of the on-line Sup- 
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plement to ref. [3]), or simply refuted as being identically 
valid irrespective of the value of an energy parameter x 
(see for instance p. 4 of the on-line Supplement to ref. 
[5]). With the help of the Euler theorem, which trans- 
forms the continued fractions into an infinite series, our 
approach turned the Schweber quantization criterion into 
an efficient computational tool. 

Our approach has been demonstrated on the examples 
of a displaced harmonic oscillator and the Rabi model. 
However, the outlined approach would work for any of the 
models of the class 1Z. For instance, the JC model and the 
single-mode spin-boson form of a generalized Rabi model 
introduced in [5], 

H Rt = cja' a + Aer 3 (a t + a) + tu 3 + fxai, (41) 

where r is a deformation parameter. Another example is 
a modified Rabi model with the interaction Hamiltonian 
Hi n t = iHX(Ti(cL ! — a). The latter arises if the dipole in- 
teraction with a Fabry-Perot cavity light mode is replaced 
by the interaction with a single-mode plane-wave field. 

We have characterized our class 7Z of models implic- 
itly in that their eigenvalue equation can be reduced 
to a difference equation ((2|) with linin^oo a n and 
lim„_ > . 00 b n /a n = 0. An interesting problem is to provide 
an explicit characterization of the models. 

Conclusions. — A general formalism has been devel- 
oped which allows to determine the spectrum of an entire 
class of quantum models as zeros of a corresponding tran- 
scendental function F(x). The function can be analyti- 
cally determined as an infinite series defined solely in terms 
of recurrence coefficients. The class of quantum models 
comprises the displaced harmonic oscillator, the Jaynes- 
Cummings (JC) model, the Rabi model, and a generalized 
Rabi model. Applications of the Rabi model range from 
quantum optics and magnetic resonance to solid state and 
molecular physics. The model plays a prominent role in 
cavity QED and circuit QED, and can be experimentally 
realized in Cooper-pair boxes, flux q-bits, in Josephson 
junctions or using trapped ions, and is of importance for 
various approaches to quantum computing. Therefore, our 
results could have implications for further theoretical and 
experimental work that explores the interaction between 
light and matter, from weak to strong interactions. The 
ease in obtaining the spectrum is of importance regarding 
recent experimental advances in preparing (ultra) strongly 
interacting quantum systems, which can no longer be reli- 
ably described by the exactly solvable JC model. The rel- 
evant computer code has been made freely available online 

* * * 
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